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Abstract 



The calibration method is used to identify some minimizers of the 
Mumford-Shah functional. The method is then extended to more general 
free discontinuity problems. 



1 Introduction 

In |^| De Giorgi introduced the name free discontinuity problems to denote a 
wide class of minimum problems for functionals of the form 

F(u):=[ f{x,u{x),\7u(x))dx + [ ip(x,u + (x),u-(x),i> u (x))dH n -\ (1) 
Ja\s u J s u 

where VL is a given bounded domain in R™ with Lipschitz boundary, /: f2 x R x 
R" — > [0, +oo] and V> : ^ x R x R x S™ -1 — > [0, +oo] are given Borel functions, 
S n_1 = {v e R™ : |u| = 1}, H 71 ^ 1 is the (n — l)-dimensional Hausdorff measure, 
and the unknown function u: Q — > R is assumed to be regular out of a (partially 
regular) singular set S u of dimension n — 1, with unit normal v u , on which u 
admits unilateral traces u + and uT . The main feature of these problems is that 
the shape and location of the discontinuity set S u are not prescribed. Thus 
minimizing F means optimizing both the function u and the singular set S u , 
which is indeed often regarded as an independent unknown. 

These problems have an increasing importance in many branches of applied 
analysis, such as image processing (Mumford-Shah functional for image seg- 
mentation) and fracture mechanics (Griffith's criterion and Barenblatt cohesive 
zone model). 

The Mumford-Shah functional was introduced in |jTof in the context of a 
variational approach to image segmentation problems (for which we refer to 
[QJ). It can be written, in dimension n, as 

F?>P(u):= [ \Vu(x)\ 2 dx + aH n - 1 (S u )+ [3 [ \u[x) - g(x)\ 2 dx , (2) 
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where g is a given function in L°°(fi) (interpreted as the grey level of the image 
to be analysed), and a > and f3 > are constants. When n = 2 (the only case 
considered in image processing), the singular set S u of a minimizer u of F^^ is 
interpreted as the set of the most relevant segmentation lines of the image. 

Using different classes of infinitesimal variations, one can show that every 
minimizer must satisfy certain equilibrium conditions, which could be globally 
called Euler-Lagrange equations for F^^ . For instance, u must satisfy the equa- 
tion Au — (3(u — g) on Q\S U , with Neumann boundary conditions on S u U <9f2. 
Moreover, there is a link between the mean curvature of S u (where defined) 
and the traces of u and Vw on the two sides of S u ; for instance, when (3 = 0, 
the mean curvature of S u must be equal to the difference of the squares of the 
norms of the traces of Vw. Additional conditions have been derived for the two- 
dimensional case. We refer the reader to jlO| and for a precise description of 
these equilibrium conditions. 

However, since F^^ is not convex, all conditions which can be derived by 
infinitesimal variations are necessary for minimality, but never sufficient. The 
purpose of this note is precisely to present a sufficient condition for minimality 
(Theorem 3.1 for F^^ and Theorem |3.4| for F), and give a few applications 



(Examples 4.1-4.8). Detailed proofs and further results will be given in the 
forthcoming paper 



2 Notation and preliminaries 

For a complete mathematical treatment of the minimum problems for the func- 
tional F considered in (|l|), we use the space SBV(Q) of special functions of 
bounded variation, introduced by De Giorgi and Ambrosio in A self-con- 
tained presentation of this space can be found in the recent book J2|, which 
contains also the complete proof of the existence of a minimizer u of Fg'P, and 
of the partial regularity of the corresponding singular set S u (the regularity of 
u on fl \ S u follows from the standard theory of elliptic equations). 

We recall that for every u S SBV(fl) the approximate upper and lower limits 
u + (x) and u~(x) at a point x G Q arc defined by 

u^x) :=±inf{iGR: lim p- n C n ({±u > t} n BJx)) = 0} , 

where B p {x) is the open ball with centre x and radius p. The singular set (or 
jump set) of u is defined by S u := {x £ Q : u~(x) < u + (x)}. It is known that S u 
is countably (TC 1 ^ 1 ,n — l)-rcctifiable and that there exists a Borel measurable 
function v u : S u — * S"^ 1 such that for H n ~ 1 -a.e. x S S u we have 

lim ^ f | U (y)- U ± (x)|dy = 0, (3) 
p-o+ p n J B ± {x) 

where B^(x) := {y e B p (x) : ±(y — x) ■ v u {x) > 0} and • denotes the scalar 
product in R™ (see 0, Theorem 4.5.9]). Condition (||) says that v u {x) points 
from the side of S u corresponding to u~(x) to the side corresponding to u + (x). 
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The gradient Du of u is a measure that can be decomposed as the sum of 
two measures Du = D a u + D s u, where D a u is absolutely continuous and D s u 
is singular with respect to the Lebesgue measure £ n . The density of D a u with 
respect to C n is denoted by Vit. Since u £ SBV(fl), for every Borel set B in SI 
we have 

{Du){B)= [ X7u(x)dx+ [ (u+(x) - u-{x)) v u (x) dH^ 1 . 
JB JBnS u 

The graph of u is defined as 

r u :— {(x,t) £ fl x K : u~(x) < t < u + (x)} . 

The characteristic function of the subgraph {(x,t) £ ft x R : t < u(x)} is de- 
noted by l u . It is defined by l u (x,t) := 1 if t < u(x), and l M (a;,t) := if 
t > u{x). It belongs to SBV(Q x R) and its gradient Dl u is a measure concen- 
trated on r u . 



3 The main results 

We fix an open subset 17 of O x R of the form 

U := {(x,t) £ Q x R : < i < t 2 (x)} , (4) 

where t± and r 2 are two continuous functions on £1 such that — oo < t±(x) < 
T2(x) < +oo for every x £ CI. 

Let F be the functional introduced in ([!]). We say that a function u £ 
SBV(n), with graph ^"-contained in U (i.e., H n (r u \ U) = 0), is a Dirichlet 
U -minimizer of -F, if < for every t> G 573^(0) with the same trace 
as u on dil and with graph 7i"-contained in U. If the inequality F(u) < F(v) 
holds for every v £ SBV(fl) with graph 7i"-contained in U, we say that u is a 
C7 -minimizer of -F. We omit [7 when U = il x R. 

The symbol will always denote a bounded Borel measurable vectorfield 
defined on U with values in R™ +1 — R™ x R, with components <j) x £ R n and 
0* £ R. The divergence of (f> is then div</>(:r,i) = div x (j) x (x, t) + dt^ix, t). 

We begin with a theorem concerning the functional Fg'@ introduced in (||). 

Theorem 3.1 Let u £ SBV(fl) with graph TL n -contained in U. Assume that 
there exists a bounded vectorfield 4> of class C 1 on U with the following proper- 
ties: 

(al) \\r(x,t)\ 2 <<j> t (x 1 t)+(3\t-g(x)\ 2 

for C n -a.e. x £ Q and for every T\{x) < t < t 2 (x); 

(a2) cjf(x,u(x))=2Wu(x) and 4>\x, u{x)) = | Vii(a;)| 2 - (3\u(x) - g(x)\ 2 
for C n -a.e. x £ Q; 
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(bl) 



(b2) 



(cl) 



(■l 2 



(j) x (x,t)dt 



< a 



for TL n 1 -a.e. x G Q and for every T\{x) < <i < i2 < T2(x); 

(x) 

4> x {x, t)dt — a v u {x) for Ti 71 -a.e. x € S u ; 

u~ (x) 

drvcj)(x,t) = for every (x,t) G U. 



Then u is a Dirichlet U-minimizer of Fg" . If, in addition, cf> x (x,t) satisfies the 
boundary condition 

(c2) lim <j) x {y,s)-v{x) =0 

(y,s)^(x,t) 

for H"^ 1 -a.e. x € dfl and for C 1 -a. e. t£ [ti(x), T2(x)], 
where v(x) is the outer unit normal to dtt, then u is a U-minimizer of Fg^ . 



A vectorfield 4> which satisfies conditions (al)-(cl) of Theorem 3.1 is called 
a calibration for the functional F^ 3 on U. If <p satisfies also (c2), it is called a 
Neumann calibration. Theorem 3.1 is an immediate consequence of the following 
lemmas. 

Lemma 3.2 Let <j) be a vectorfield which satisfies conditions (al) and (bl) of 
Theorem 3.1. Then for every u S SBV(Q) with graph TV 1 -contained in U we 
have 

Ff\u)> f <p-d(Dl u ) . (5) 
Ju 

Moreover, equality holds in ^) for a given u if and only if conditions (a2) and 
(b2) of Theorem 3.1 are satisfied. 

The next lemma is a consequence of the divergence theorem. 

Lemma 3.3 Suppose that <f> is of class C 1 and that div0 = on U . Then 

dp-d{Dl u ) = [ <P-d(Dl v ) (6) 



for every pair of functions u, v in BV(Sl) with the same trace on <9f2 and with 
graphs TL n -contained in U . If, in addition, </> satisfies condition (c2) of Theo- 
rem 3.1, then (&) holds for every pair of functions u, v in BV(ft) with graphs 



H. n -contained in U. 

As a matter of fact, the method of calibrations can be easily adapted to the 
functional F defined in ([l]). 

Theorem 3.4 Let u G SBV(fl) with graph HI 1 -contained in U . Assume that 
there exists a bounded vectorfield if) of class C 1 on U with the following proper- 
ties: 
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(al) 4> x {x,t) -v < 4>\x,t) +f(x,t,v) forC n -a.e. xeto, 
for every T\[x) < t < T2(x), and for every v G R"; 

(a2) (j?{x, u{x)) ■ Vu(x) = ^(x, u(x)) + f(x, u(x),Vu(x)) 
for C n -a.e. x € to; 

(bl) v \ (j) x (x,t)dt<ip(x,ti,t 2 ,v) for H^-a.e. x e to, 
Jti 

for every Ti(x) < t\ < t 2 < T2(x), and for every v 6 S" _1 ; 

(b2) v u {x)- / 4> x (x,t)dt = ip(x, u~{x),u + (x),v u (x)) 

forU^-a.e. x £ S u ; 
(cl) dW(p(x,t) = for every (x,t) € C/. 

T/ien u is a Dirichlet U-minimizer of F. If (p x (x,t) satisfies also the boundary 



condition (c2) of Theorem 3.1, then u is a U-minimizer of F. 



Remark 3.5 We note that in Theorem 3.4 there is no regularity or convexity 
hypothesis on / or ip. If f*(x, t, v*) is the the convex conjugate of fix, t, v) with 
respect to v, condition (al) is equivalent to 

(al') f*{x,t,cjf{x,t))<^{x,t) 

for £ n -a.e. x € to and for every t\ (x) < t < t 2 (x) . 

If this condition is satisfied, and f(x,t,v) is convex and differentiable with 
respect to v, then condition (a2) is equivalent to 



, d> x (x, u(x)) = d v f (x, u(x), Vu(x)) 
(a2') ^ A W y " for £"-a.e. x € to. 

</>*(£, u(aO) = f*(x,u(x),4> x (x,u(x))) 



Remark 3.6 In Theorems 3.1 and 3.4 the hypothesis that <f> is of class C 1 is 



too strong for many applications. It is used only in Lemma 3.3 and it can be 
relaxed in several ways (see [Q for details). For instance, one may consider 
piecewise C 1 vectorfields, which may be discontinuous along sufficiently regular 
interfaces. In this case the divergence-free condition (cl) must be understood in 
the distributional sense, i.e., the pointwise divergence vanishes (where defined) 
and the normal component of <f> is continuous across the discontinuity surfaces. 
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4 Some examples 



The following examples show that the calibration method is very flexible, and 
can be used to prove the minimality of a given function u in many different situ- 
ations. In the first examples we will consider only the "homogeneous" functional 
F a :— Fg-°, in which the lower order term (3 J n \u — g\ 2 dx vanishes. 

Example 4.1 (Affine function in one dimension) Let n := 1, fi := ]0, a[, 

and u(x) := Ax, with A > 0. It is easy to see that u is a Dirichlet minimizer of 
F a if and only if a\ 2 < a. In this case a calibration is given by the piecewise 
constant function 

J(2A,A 2 ), if $x<t<${x + a), 
<j>{x,t) ■= < (7) 
[ (0, 0), otherwise. 

Another calibration is given by 

f (2|,(|) 2 ), ifO<t<Xx, 

4>(x, t) := | (2^=*, (^) 2 ) , if Ax < t < Xa, (8) 

[ (0, 0), otherwise. 

If aX 2 > a, then the function u(x) :— Xx is not a Dirichlet minimizer of F a , 
but it is still a Dirichlet U -minimizer with 

U := {(x,t) e}0 7 a[x R, : Xx - < t < Xx + ^} . 

A calibration on U is given by <j>(x, t) :— (2A, A 2 ). 



Example 4.2 (Jump in one dimension) Let n := 1, f2 := ]0, a[, u(x) := 
for < x < c, and u{x) := h for c < x < a, with < c < a and h > 0. It is easy 
to see that u is a Dirichlet minimizer of F a if and only if act < h 2 . In this case 
two different calibrations are given by (@) and (§) with A 

Suppose now that aa > h 2 . Let e > be a constant such that 2e J ^^/2ae < h, 

let 

!— £, if x < c, 

-e+|(a;-c), ifc<x<c + e, 
h — s, if c + e < x, 

let T2(x) = Ti(x + e) + 2s, and let U be the open set defined by (^). Then u is 
a Dirichlet [/-minimizer of F a , and a calibration on U is given by the piecewise 
constant function 

!(2A, A 2 ), ifc-e<x<c + e and 
s + |(x - c + e) < t < £ + |(x - c + e) + 
(0,0), otherwise, 

where A > is any constant such that e + eX + < h — e, for instance 
A = ,/a/V2e. 
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Example 4.3 (Harmonic function) Let fi be a bounded domain in R", n 
arbitrary, and let u be a harmonic function on f2. As pointed out by Chambollc 
[|J, u is a Dirichlet minimizer of F a if 

oscusup|Vw| < a, (9) 

where oscq u := supn u — inf^ u. Note that for n = 1 this condition reduces to 
the constraint aX 2 < a of Example 4.1. Inspired by the one dimensional case 
(see (0)), we construct the calibration 

/ (2V U (x),|Vu(x)| 2 ), if §( U (x) + m) < t < \{u{x) + M), 
4>{x,t) := < - (10) 

[ (0, 0), otherwise, 

where m :— irrfn u and M :— supo w. Another calibration (see (^)) is given by 
f (2^^^), (-^-)2|Vu( a ;)| 2 ), if m < t < «(*), 
^*):=| ( 2 l^)V U ( a; ),( 1 ^ y ) 2 |Vu( a; )| 2 ), if^)<t<M, 
[ (0, 0), otherwise. 
If (|9|) is not satisfied, u is still is a Dirichlet [/-minimizer of F a , for 
£7 := {(x,t) e n x R : w(x) - f {Vuix)^ 1 < t < u{x) + f ^^(a;)!- 1 } , (11) 
and a calibration in [/ is given by 4>(x,t) :— (2Vu(x), |Vm(x)| 2 ). 



Example 4.4 (Pure jump) Let n > 2 and let f2 := ]0, a[ x 1/, where 1/ is a 
bounded domain in R ,l_1 with Lipschitz boundary. Denoting the first coordi- 
nate of a: by x\, let u(x) :— for < x\ < c, and u(x ) :— /i for c < xi < a, with 
< c < a and ft, > 0. Using the results of Example [b^ it is easy to see that u 
is a Dirichlet minimizer of F a if aa < /i 2 . In this case two different calibrations 
can be constructed in the following way: the projection of these calibrations 
onto the (xi, i)-plane are given by (Q) and (^), with A = y/a/y/a and x replaced 
by xi, while all other components of these calibrations vanish. 

If aa > h 2 , it may happen that u is still a Dirichlet minimizer of F a . For 
instance, if n = 2 and V = ]0,b[, with bair < 2h 2 , a different calibration has 
been constucted in 0]. Therefore u is a Dirichlet minimizer of F a even if aa is 
very large with respect to h 2 , provided that ba is small enough. 

Arguing as in the last part of Example [4.2| one can prove that for every a 
and V there exists an open set U of the form (14), containing r u , such that u is 
a Dirichlet [/-minimizer of F a . 



Example 4.5 (Triple junction) Let n := 2, let £1 := B(0, r) be the open ball 
with radius r > centered at the origin, and let u be given, in polar coordinates, 
by u(p,9) := a for < 9 < |7r, u(p,6) := b for |-7r < 9 < |7r, and u(p,9) := c 
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for |-7r < 9 < 2tt, where a, b, and c are distinct constants. Thus S u is given by 
three line segments meeting at the origin with equal angles. If 

2ar < min{|a - b\ 2 , \b - c| 2 , |c - a| 2 } , (12) 

then u is a Dirichlet minimizer of F a . To construct a calibration, it is not 
restrictive to assume a < b — < c. Inspired by the one dimensional case 



described in Example 4.2, we take e± := (±V3/2, -1/2), and A > such that 
4f + j < min{— a, c} (which is possible by (|l2|)), and we define the calibration 

by 

(Ae+, A 2 /4), if |(r + x ■ e+) < t < |(r + x ■ e+) + f , 
{x, t) ■= \ (Ae_, A 2 /4), if A(_ r + a, . e _) _ s < t < A ( _ r + x . e _ ); (13 ) 

(0,0), otherwise. 

If ar is much larger than min{|a — b\ 2 , \b — c| 2 , |c — a| 2 }, it is easy to construct 
a comparison function v with the same boundary values as u and such that 
F a (v) < F a {u). This shows that in this case u is not a Dirichlet minimizer. 

However, for every value of the parameters a, r, a, 6, c, one can construct a 
suitable neighbourhood {J of the graph r u , of the form (|j), such that a variant 
of (|i~3] ) is a calibration in U, and therefore u is a Dirichlet [/-minimizer of F a . 
We refer to for the details. 

We consider now the functional F^^ , with /? > 0. 

Example 4.6 (Solution of the Neumann problem) Let be a bounded 
open set in R" with boundary of class C 1,£ for some e > 0, and let u be the 
solution of the Neumann problem 

{Am = 8(u — a) on fi, 
(14) 
= o on on, V ; 

is 



with P > and g € L°°(f2). Assume that condition (g|) of Example 
satisfied. Then u is a minimizer of F^*'' 3 . If the strict inequality holds in (|), 
then w is the unique minimizer. A Neumann calibration <f>(x, t) is given by 



2 • 



( ,/3| f _ _ ^ f + uM _ ff(x) | 2)i if t - «M < 

(2V«(a;) ) IV^)! 2 - /3|t - sWI 2 + P\t ~ u(x)f), if f < t - *f < f , 

( ,/?|f -HM|2_£|M + «M_ g(a; )|2) j if M <f _ «M 

where m := infn u and M := supo u. 

If @ is not satisfied, u is still is a [/-minimizer of -f 7 '^'' 3 , where U is the open 
set defined by (fTO). A Neumann calibration on U is given by 



<j>{x,t) := (2Vw(a;),|VM(a;)| 2 - /3|t - ,g(a;)| 2 + (3\t - u(x)\ 2 ) . 
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The hypothesis that dft is of class C 1,6 is used only to obtain the boundary 
condition (c2) of Theorem 3.1, which, in this case, becomes 



lim Vu(y) ■ v(x) = for TT^-a.e. x £ dft . (15) 

It is clear that ( |i~5| ) is still true if for H™ _1 -a.e. x G dft there exists an open 
neighbourhood V x of x in R™ such that V x n <9fi is a manifold of class C 1,e (see 
[|[ Theorem 7.5.2]). Therefore the result of this example is true also when ft is 
polyhedral. 

In the next examples we construct a calibration for Fg'" when the parameter 
(3 is large enough. 

Example 4.7 (Smooth g and large (3) Let ft be a bounded open set in R" 
with smooth boundary, and let g S C 2 (f2). There exists a constant /3o > 0, 
depending on g and a, such that for every f3 > (3q the solution u of the Neumann 

is the unique minimizer of F^ 13 . A Neumann 
calibration is constructed in Jl| . 

This shows that the minimizer of Fg^ is smooth, provided that g is smooth 
and (3 is large enough. Therefore the solution of the image segmentation problem 
(n = 2) based on the minimization of Fg^ has an empty set of segmentation 
lines if the "grey level" function g is smooth and the parameter f3 in the fidelity 
term (3 J Q \u — g\ 2 dx is large. 

Example 4.8 (Function g with only two values) Let ft be an open set in 
R™ and let E be a compact set contained in ft with boundary of class C 2 . Let 
g{x) := a for x € E and g{x) :— b for x € Q \ E, with a b. There exists 
a constant /3q > 0, depending on g and a, such that for every (3 > /3q the 
function u :— g is the unique minimizer of Fg'@. To construct a calibration, it 
is not restrictive to assume a < b. We take a C 1 vectorfield v. £1 — > R™ with 
compact support in ft such that \v(x)\ < 1 for every x S ft and v(x) is the 
outer unit normal to dE for every x G Then we set (f> x (x,t) — <j(t)v(x), 
where a is a fixed positive smooth function with integral equal to a and support 
contained in ]a, b[. We see that conditions (bl), (b2), and (c2) of Theorem |3.1| 
are satisfied by construction. It remains to choose (/>* so that (al), (a2), and 
(cl) hold. Condition (a2) forces us to set (^{x^t) — for t — g{x), while (cl) 
gives d t (j) t {x 1 t) = —a(t)div x v(x). These two conditions determine 0*(a;,i) at 
every point (x, t). It is then easy to see that (al) holds if (3 is large enough. We 
refer to [Q for the details. 

This example shows that, if g £ SBV(ft) has only two values, and S g is 
smooth enough, then the minimizer of the Mumford-Shah functional Fg'& re- 
constructs g exactly, when (3 is large enough. 

Recently the following question has been studied by using the calibration 
method: is it true that a function u is a (Dirichlet) minimizer of Fg'P, if it 
satisfies the Euler-Lagrange equations and the domain fl is sufficiently small? 



problem (p4) of Example 4.6 
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For the moment we have only a partial answer. In 0] we have considered the 
case where n := 2 and S u is a line segment joining two points of the boundary of 
f2. If u satisfies the Euler-Lagrange equations for the "homogeneous functional" 
F a := Fg'°, then for every xq £ S u there exists an open neighbourhood f2o 
of Xq, contained in f2, such that u is a Dirichlet minimizer of F a in Qq. The 
minimality is proved by constructing a complicated calibration on Slo x 

This result has been extended in j§] to the case where S^i is an analytic curve 
joining two points of dft. The (more difficult) construction of the calibration 
presented in this paper shows that one can take the same set fio for every 
xq G S u ; in other words, one can take as f2o a suitable tubular neighbourhood 
of S u - Moreover, it is proved in B that an additional condition on u and S u 
implies that u is a Dirichlet [/-minimizer for a suitable open neighbourhood U 
of the graph r u . A counterexample (where S u is a line segment joining two 
points of dtt) shows that this is not always true when u is just a solution of 
the Euler-Lagrange equations with S u ^ 0, in contrast to the case S u = (see 
Example 
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